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We developed a three-center phenomenological model, able to explain qualitatively 
the recently obtained experimental results concerning the quasimolecular stage of 
a light-particle accompanied fission process. It was derived from the liquid drop 
model under the assumption that the aligned configuration, with the emitted par- 
ticle between the light and heavy fragment, is reached by increasing continuously 
the separation distance, while the radii of the heavy fragment and of the light par- 
ticle are kept constant. In such a way, a new minimum of a short-lived molecular 
state appears in the deformation energy at a separation distance very close to the 
touching point. This minimum allows the existence of a short-lived quasi-molecular 
state, decaying into the three final fragments. The influence of the shell effects is 
discussed. The half-lives of some quasimolecular states which could be formed in 
the 10 Be and 12 C accompanied fission of 252 Cf are roughly estimated to be the 
order of 1 ns, and 1 ms, respectively. 

1 Introduction 

The light particle accompanied fission was discovered in 1946, when the 
track of a long-range particle (identified by Farwell et al. to be 4 He) almost 
perpendicular to the short tracks of heavy and light fragments was observed in 
a photographic plate. The fission was induced by bombarding 235 U with slow 
neutrons from a Be target at a cyclotron. The largest yield in such a rare pro- 
cess (less than one event per 500 binary splittings) is measured for a-particles, 
but many other light Jiuclei (from protons to oxygen or even calcium isotopes) 
have been identified □ in both induced- and spontaneous fission phenomena. 
If Ai and Ai are the mass numbers of the heavy fragments (assume A\ > A2), 
then usually the mass of the light particle A3 << A%. The "true" ternary 
fission, in which A\ ~ A2 ~ A3, has not yet been experimentally detected. 
Many properties of the binary fission process have been explained a within 
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Figure 1. Gamma-ray energy spectrum of 10 Be accompanying fission of 252 Cf (top) and the 
simulated (with the code GEANT) pulse-height distributions for a (not-Doppler-broadened) 
3.37 MeV 7 emission from 10 Be at rest (middle) as well as for 10 Be in flight (bottom) from 
the source to the detector (duration of about 1 ns). The flight-path of the order of 10 cm. 
One uses Nal scintillator detectors. Data from Fig. 6 of Ref. 13. 



the liquid drop model (LDM); others like the asymmetric mass distribution 
of fragments and the ground state deformations of many nuclei, could be 
understood only after adding the contribution of shell effects. Qi3 As it was 
repeatedly stressed (see H and the references therein), shell effects proved also 
to be of vital importance for cluster radioactivities predicted in 1980. 

The total kinetic energy (TKE) of the fragments, in the most frequently 
detected binary or ternary fission mechanism, is smaller than the released 
energy (Q) by about 25-35 MeV, which is used to produce deformed and 
excited fragments. These then emitt neutrons (each with a binding energy 
of about 6 MeV) and 7-rays. From time to time a "cold" fission mechanism 
is detected, in which the TKE exhausts the Q-value, hence no neutrons are 
emitted, and the fragments are produced in or near their ground-state. The 
first experimental evidence for cold binary fission in which its TKE exhaust Q 
was reported E in 1981. Larger yields were measured □ in trans-Fm [Z > 100) 
isotopes, where the phenomenon was called bimodal fission. 

The correlated fragment pairs in cold ternary (a- and 10 Be accompanied 
spontaneous fission of 252 Cf) processes were only recently discovered, E30 
by measuring triple 7 coincidences in a modern large array of 7-ray detec- 
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Figure 2. The cumulated 7-ray energy spectrum of 10 Be accompanying fission of 252 Cf 
(top) and the corresponding background (bottom). The spectrum of the 3.362 7 is not- 
Doppler-broadened, suggesting an emission from 10 Be at rest. The stopping time of 1(> Be 
in the absorber mounted around the source is of the order of 1 ps. Data from Fig. 2 of Ref. 
11. Coincidence spectra obtained with Ge detectors in GAMMASPHERE. 



tors (GAMMASPHERE). The fragments are identified by their 7-ray spectra. 
Among other new asnects of the fission process seen for the first time with 
this new technique, E3'E3 one should mention the double fine structure, and 
the triple fine structure in binary and ternary fis sion 

A particularly interesting feature, observed ElE3 both in 10 Bc- and 12 C 
accompanied cold fission of 252 Cf is related to the width of the light particle 
7-ray spectrum (see Figs, [l] and |J). For example, the 3.368 MeV 7 line of 
10 Be, with a lifetime of 125 fs is not Doppler-broadened, as it should be if 
it would be emitted when 10 Be is in flight (taking about 1 ns to reach the 
detector). A plausible suggestion was made, that the absence of Doppler 
broadening is related— to a trapping of 10 Be in a potential well of nuclear 
molecular character. Eil 

Quasi-molecular configurations of two nuclei have been suggested as a 
natural explanation for the resonances measured E-3 in 12 C+ 12 C scattering 
and reactions. There are also other kinds of such binary molecules (see lJ and 
references therein), like spontaneously fissioning shape- isomers. The above 
mentioned experiments can be considered as the first evidence for a more 
complex quasi-molecular configuration of three nuclei. The purpose of the 
present lecture is to show, within a phenomenological three-center model, 
that a minimum which could explain the existence of these quasi-molecules 
is produced in the potential barrier, when the formation of the light particle 
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occurs in the neck between the two heavier fragments. In this way we extend 
to ternary fission our unified approach of cold fission, cluster radioactivities, 
and a-decay. □ 



2 Shape Parametrization 

The shape parametrization with cp£ deformation parameter as follows has 
been suggested from the analysis Ii3 of different aligned and compact con- 
figurations of fragments in touch. A lower potential barrier for the aligned 
cylindrically-symmetric shapes with the light particle between the two heavy 
fragments, is a clear indication that during the deformation from an initial par- 
ent nucleus to three final nuclei, one should arrive at such a scission point. In 
order to reach this stage we shall increase continuously the separation distance, 
i?, between the heavy fragments, while the radii of the heavy fragment and 
of the light particle are kept constant, R± = constant, R3 = constant. Unlike 
in the previous work, we now adopt the following convention: A\ > Ai > A3. 
The hadron numbers are conserved: A\ + A2 + A3 = A. 

At the beginning (the neck radius p n eck > R3) one has a two-center 
evolution (see Fig. ||) until the neck between the fragments becomes equal 
to the radius of the emitted particle, p nec k — p(z s i) \r^r ov3 = R3- This Eq. 
defines R OV 3 as the separation distance at which the neck radius is equal to R3. 
By placing the origin in the center of the large sphere, the surface equation 
in cylindrical coordinates is given by: 

2 / R\-z 2 , -Ri < z < z s i , , 

Ps \R\-{z-Rf , z sl <z<R + R 2 1 ' 



• O i • 

•3 » 



Figure 3. Evolution of nuclear shapes during the deformation process from one parent 
nucleus 252 Cf to three separated fragments 146 Ba, 10 Be, and 96 Sr. In the upper part 
the binary stage is illustrated; the separation distance increases from Ri to R OV 3, passing 
through R m i n ib an d R m in2b values. In the middle, the ternary stage of the process develops 
by forming the third particle in the neck. The quasi-molecular shape, at which R = R m i n — t 
is the intermediate one in this row. At the bottom the fragments are separated. 
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Then for R > R OV 3 the three center starts developing by decreasing pro- 
gressively with the same amount the two tip distances h\ + /131 = /132 + h 2 . 
Besides this constraint, one has as in the binary stage, volume conservation 
and matching conditions. The R 2 and the other geometrical quantities arc 
determined by solving numerically the corresponding system of algebraic equa- 

1 /3 

tions. By assuming spherical nuclei, the radii are given by Rj = 1.2249A •' fm 

(j = 0,1,3), R 2f = 1.2249A2 73 with a radius constant r = 1.2249 fm, from 
Myers-Swiatecki's variant of LDM. Now the surface equation can be written 
as 



and the corresponding shape has two necks and two separating planes. Some 
of the important values of the deformation parameter R are the initial distance 
Ri = Rq — Ri, and the touching-point one, R t = R\ + 2i? 3 + R 2 j. There is 
also R OV 3, defined above, which allows one to distinguish between the binary 
and ternary stage. 

3 Deformation Energy 

According to the LDM, by requesting zero energy for a spherical shape, the 
deformation energy, E U (R) — E°, is expressed as a sum of the surface and 
Coulomb terms 



where the exponent u stands for uniform (fragments with the same charge 
density as the parent nucleus), and refers to the initial spherical parent. In 
order to simplify the calculations, we initially assume the same charge density 
Pie = P2e = P3e = Poe, and at the end we add the corresponding corrections. 
In this way we perform one numerical quadrature instead of six. For a spher- 
ical shape E° s = a s (l - kI 2 )A 2 / 3 ; I = (N - Z)/A; E° c = a^A' 1 / 3 , where 
the numerical constants of the LDM are: a s = 17.9439 MeV, k = 1.7826, 
a c = 3e 2 /(5r ), e 2 = 1.44 MeV-fm. 

The shape-dependent, dimensionless surface term is proportional to the 
surface area: 




-Ri < z < z sl 

Zsl < Z < Z s2 

z s2 < z<R + R 2 



(2) 



E u def {R) = E° S [B S (R) 1] + E° C [B C (R) 1] 



(3) 
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Figure 4. The liquid drop model deformation energy versus separation distance for the 
20 O accompanied cold fission of 252 Cf with 132 Sn and 100 Zr heavy fragments. In order 
to simplify the numerical calculations we start by assuming the same charge density of 
the fragments. One can see the effect of successive corrections taking into account the 
experimental Q-value and the difference in charge density. Similar curves for the binary 
fission posses a narrower fission barrier. The new minimum appears in the shaded area 
from R v3 to Rt. 



where y = y[x) is the surface equation in cylindrical coordinates with -1, +1 
intercepts on the symmetry axis, and d = (z" — z')/2i?o_is the seminuclear 
length in units of Rq. Similarly, for the Coulomb energy L3 one has 



5d 5 f f 
B c = — — dx dx'F(x, x') 

o7T J J 



(5) 



F(x,x') = {y yi [(K-2D)/3}- 

2^ + yl)-(x-x'f + \{x-x') (gf 
2 x — x' dy 2 



K\y 2 yl/Z 



dyf djf 
dx 

x — x' dy 2 



dx' 



W (6) 



2 dx 

K, K' are the complete elliptic integrals of the 1st and 2nd kind 

tt/2 tt/2 

K{k)= f (1 - k 2 sm 2 t)-^ 2 dt; K'(k)= [ (1 - k 2 smh) x ' 2 dt (7) 



and a 2 p = {y + Vl f + (x - x 1 ) 2 , k 2 = 4 yyi /a 2 p , D = (K - K')/k 2 . 

The new minimum, which can be seen in Fig. |4| at a separation distance 
R = Rmin-t > Rov3, is the result of a competition between the Coulomb- 
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and surface energies. At the beginning (R < R m i n ~t) the Coulomb term is 
stronger, leading to a decrease in energy, but later on (R > i? m i„_ t ) the light 
particle formed in the neck posses a surface area increasing rapidly, so there 
is also an increase in energy up to R — R t ■ 

Now let us analyse the influence of various corrections, which could in 
principle alter this image. After performing numerically the integrations, we 
add the following corrections: for the difference in charge densities reproducing 
the touching point values; for experimental masses reproducing the Q exp - 
value at R = Ri, when the origin of energy corresponds to infinite separation 
distances between fragments, and the phenomenological shell corrections SE 

E LD (R) = E% ef (R) + (Q th - Q exp )f c (R) (8) 

where f c (R) = {R - R t )/{R t - R t ), and 

3 3 

Q th = E° s +E° C - J^iEsi + E° Ci ) +6E°-J2 ^ (9) 
i i 

The correction increases gradually (see Fig. ^ and Fig. ||) with R up to Rt and 
then remains constant for R > R t . The barrier height increases if Q eX p ^ Qth 
and decreases if Q exp > Qth- In this way, when one, two, or all final nuclei have 
magic numbers of nucleons, Q exp is large and the fission barrier has a lower 
height, leading to an increased yield. In a binary decay mode like cluster 
radioactivity and cold fission, this condition is fulfilled when the daughter 
nucleus is 208 Pb and 132 Sn, respectively. 

4 Shell Corrections and Half-lives 

Finally we also add the shell terms 

E(R) = E LD {R) + SE(R) - SE° (10) 

Presently there is not available any microscopic three-center shell model 
reliably working for a long range of mass asymmetries. This is why we use 
a phenomenological model, instead of the Strutinsky's method, to calculate 
the shell corrections. The model is adapted after Myers and Swiatecki. u 
At a given R, we calculate the volumes of fragments and the corresponding 
numbers of nucleons Zi(R), Ni(R) {i — 1,2,3), proportional to the volume of 
each fragment. Fig. || illustrates the evolution of shapes and of the fragment 
volumes. Then we add for each fragment the contribution of protons and 
neutrons 

6E(R) =J2SE t (R) = Y^[SEpi(R) + SE m (R)} (II) 

i i 
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Figure 5. The liquid drop model, Eld, the shell correction, 8E, and the total deformation 
energies, E, for the 10 Be accompanied cold fission of 252 Cf with 146 Ba and 96 Sr heavy 
fragments. The new minimum appears in the shaded area from RovS to Rt. 



which are given by 



where 



SEpi = Cs(Zi); 8E m = Cs{N i ) 



s{Z) = F(Z)/[(Z)~ 2 / 3 ] - cZ 1 / 3 



F(n) = - 



N 



5/3 



N 



5/3 



Ni - Ni. 



,5/3 . 



n: 



5/3 



(12) 

(13) 
(14) 



in which n G (JVj_i,iVi) is either a current Z or N number and A^_i,iVj 
are the closest magic numbers. The constants c = 0.2, C = 6.2 MeV were 
determined by fit to the experimental masses and deformations. The variation 
with R is calculated as 



SE(R) = -\Y^{s(N i ) + s(Z i : 



Ri 



(15) 



where Li{R) are the lengths of the fragments along the axis of symmetry, at 
a given separation distance R. During the deformation, the variation of sep- 
aration distance between centers, R, induces the variation of the geometrical 
quantities and of the corresponding nucleon numbers. Each time a proton or 
neutron number reaches a magic value, the correction energy passes through 
a minimum, and it has a maximum at midshell (see Fig. | and Fig. g). The 
first narrow minimum appearing in the shell correction energy SE in Fig. ||, 
at R — Rminib — 2.6 fm, is the result of almost simultaneously reaching the 
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Figure 6. The shell correction energy variation with the separation distance for two exam- 
ples of 10 Be-, two of 12 C accompanied cold fission, compared to the "true" ternary fission 
(in nearly three identical fragments) of 252 Cf. The three partners are given. The two 
vertical bars on each plot show the positions of RovZ and of Rt. 

magic numbers Z\ = 20, Ni = 28, and Z2 = 82, N2 — 126. The second, more 
shallower one around R = R m in2b — 7.2 fm corresponds to a larger range of R- 
values for which Z\ = 50, N\ = 82, Z2 = 50, N2 = 82 are not obtained in the 
same time. In the region of the new minimum, R = R m i n ^ t , for light-particle 
accompanied fission, the variation of the shell correction energy is very small, 
hence it has no major consequence. One can say that the quasimolecular min- 
imum is related to the collective properties (liquid-drop like behavior). On 
the other side, for "true" ternary process (see the bottom part of Fig. |J) both 
minima appear in this range of values, but no such LDM effect was found 
there. In order to compute the half-life of the quasi-molecular state, we have 
first to search for the minimum E mm in the quasimolecular well, from R ov ^ 
to Rt , and then to add a zero point vibration energy, E v : E qs = E m i n + E v . 

The half-life, T, is expressed in terms of the barrier penetrability, P, 
which is calculated from an action integral, K, given by the quasi-classical 
WKB approximation 



T= ^eB ] P = exp{-K) (16) 
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Table 1. Calculated half-lives of some quasi- molecular states formed during the ternary 
fission of 252 Cf. 



Particle 


Fragments 


Qexp 

(MeV) 


K 


logT(s) 


10 Be 


132 Sn 


110 Ru 


220.183 


19.96 


-11.17 




138 Te 


104 Mo 


209.682 


25.23 


-8.89 




138 Xe 


104 Zr 


209.882 


26.04 


-8.54 




146 Ba 


96 Sr 


201.486 


22.98 


-9.86 


12 C 


147 La 


93 Br 


196.268 


39.80 


-2.56 




142 Ba 


98 Kr 


199.896 


42.71 


-1.30 




140 Tc 


ioo Zr 


209.728 


38.21 


-3.25 




132 Sn 


108 Mo 


223.839 


31.46 


-6.18 



where h is the Planck constant, and 




R a 



in which i? a ,i?f, are the turning points, defined by E(R a ) = E(Rb) = E qs 
and the nuclear inertia is roughly approximated by the reduced mass fi = 
m[(AiA2 + A$A)/{A\ + A 2 )], where m is the nucleon mass, log[(Mn2)/2] = 

-20.8436, loge = 0.43429 and ^8m/h 2 = 0.4392 MeV-^xfm" 1 . 

The results of our estimations for the half-lives of some quasimolecular 
states formed in the 10 Be- and 12 C accompanied fission of 252 Cf are given in 
Table 1. They are of the order of 1 ns and 1 ms, respectively, if we ignore the 
results for a division with heavy fragment 132 Sn, which was not measured due 
to very high first excited state. Consequently the new minimum we found can 
qualitatively explain the quasimolecular nature of the narrow line of the 10 Be 
7-rays. 

It is interesting to note that the trend toward a split into two, three, 
or four nuclei (the lighter ones formed in a long neck hetween the heavier 
fragments) has been theoretically demonstrated by Hill, 113 who investigated 
the classical dynamics of an incompressible, irrotational, uniformly charged 
liquid drop. No mass asymmetry was evidenced since any shell effect was 
ignored. 

In conclusion, we should stress that a quasimolecular stage of a light- 
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particle accompanied fission process, for a limited range of sizes of the three 
partners, can be qualitatively explained within the liquid drop model. 
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